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INTRODUCTION 

The  aim  of  this  ongoing  project  is  to  implement  the  new  speckle 
imaging  technique  proposed  in  "Speckle  imaging,  photon  by  photon" 

(L.  Mertz,  1979,  Applied  Optics .  18 ,  611-614;  enclosed  as  Appendix  A). 
Successful  speckle  imaging  would  permit  much  higher  resolution  images 
of  the  sun  than  are  presently  available.  Such  images  could  then  be 
used  for  detailed  interpretation  of  the  solar  magnetic  field  structure. 

Our  proposed  speckle  imaging  technique  is  contingent  on  observational 
data  from  a  2-D  photon  counter  that  preserves  the  temporal  order  of  the 
incident  photons,  as  opposed  to  the  images  that  are  provided  by  customary 
observing  devices. 

PRIOR  2-D  PHOTON  COUNTERS 

For  a  perspective,  the  two  articles,  "Low  Light  Level  Detectors  for 
Astronomy"  (P.  B.  Boyce,  1977,  Science,  198 ,  145-148^  and  "Detectors 
for  the  Space  Telescope"  (T.  Kelsall,  1978,  SPIE  Proc. ,  i  .3 ,  27-35) 
review  the  state-of-the-art  concerning  photon  counting  image  detectors. 
Those  articles  describe  concern  about  sensitivity  to  individual  photons 
as  discrete  events  and  there  is  considerable  emphasis  on  CCD  devices, 
but  negligible  consideration  for  preserving  the  temporal  ordering  of  the 
photons.  Two  of  the  devices  discussed  in  the  latter  article,  the  PHOTOCON 
and  the  MAMA  detectors,  indeed,  preserve  that  order.  However,  neither  of 
those  detectors  is  commercially  available  or  even  functioning  on  a  routine 
observational  basis,  perhaps,  because  both  include  the  position  discrimi¬ 
nation  as  a  special  internal  part  of  a  microchannel-plate  image  intensi- 
fier  tube.  The  internal  part  senses  the  position  of  the  electron  bursts 
emanating  as  events  from  the  microchannel  plate.  The  internal  location 
of  the  position  sensor  requires  that  the  construction  effort  be  delegated 
to  a  vendor  that  is  both  capable  and  willing  to  fabricate  the  special  tube. 
It  also  means  that  the  device  must  function  completely  satisfactorily  at 
the  outset,  because  adjustments  or  repairs  require  breaking  the  vacuum 
seal  and  dismantling  the  tube. 

Rather  than  hope  for  the  availability  of  either  of  those  detectors, 
a  more  expedient  and  more  cost-effective  course  seemed  to  construct  a 


detector  by  providing  external  position  sensing  of  the  photon  bursts 
emanating  as  events  from  the  phosphor  screen  of  an  almost  standard 
commercial  microchannel  plate  image  intensifier  tube. 

ANALOG  OPTICAL  POSITION  SENSING  DEVICE 

For  further  perspectives  relating  to  the  eventual  instrument,  the 
simplest  arrangement  would  involve  three  photomultipliers  viewing  the 
phosphor  screen  of  an  image  intensifier  as  shown  in  Figure  1  (where  only 
two  of  the  three  photomultipliers  are  shown) .  The  vignetting  screen 
partially  shades  one  photomultiplier  so  that  the  pulse  height  ratio  of 
the  shaded  to  the  unshaded  photomultiplier  is  a  measure  of  the  positional 
coordinate  of  the  event  on  the  phosphor  screen.  The  vignetting  screen 
for  the  third  photomultiplier  (not  shown)  is  oriented  in  the  other  direc¬ 
tion  to  furnish  the  second  coordinate  of  the  position.  Ratiometric 
analog-to-digital  converters  would  be  suitable  for  those  signal  measure¬ 
ments.  The  main  difficulty  of  this  basically  analog  scheme  is  that  1% 

U 

precision  of  position  measurement  statistically  requires  at  least  10 
photons  registered  on  each  photomultiplier.  Considering  the  quantum 
efficiency  and  the  collection  efficiency  of  the  system,  the  image  intensi¬ 
fier  tube  probably  would  need  a  gain  of  at  least  10^,  which  is  really 
high.  Other  difficulties  are  that  good  electronic  stability  is  needed 
and  that  fast  A-D  converters  are  expensive. 

DJJ :  1TAL  OPTICAL  POSITION  SENSING  DEVICE 

Those  difficulties  led  to  a  completely  digital  approach,  the  basic 
idea  being  due  to  Dr.  C.  Papaliolios.  The  idea  is  that  multiple  images 
of  the  phosphor  screen  should  be  formed  with  suitable  optics.  A  bar  mask 
corresponding  to  one  bit  of  the  event  position  is  placed  over  each  of  the 
images  along  with  a  photomultiplier.  A  tvpical  set  of  masks  is  shown  in 
Figure  2.  An  extra  photomultiplier  with  a  totally  open  mask  is  included 
so  that  for  a  2m  x  2°  pixel  picture,  m+n-t-l  images  with  masks  and  photo¬ 
multipliers  are  needed.  The  open  photomultiplier  senses  the  occurrence 
of  a  photon  burst  event;  coincident  pulses  at  any  of  the  other  photo¬ 
multipliers  indicate  that  their  corresponding  bits  should  be  ones. 
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absence  of  coincident  pulses  indicate  that  their  corresponding  bits 
should  be  zeroes.  Thus  the  positions  are  immediately  in  digital  form. 
Incidently,  the  masks  need  not  be  cartesian  or  even  rectangular  format. 
For  example,  they  may  be  arranged  to  conform  with  the  s  shaped  orders 
of  an  echelle  spectrograph  format. 

A  couple  of  initial  criticisms  of  the  idea  were  that  no  centroiding 
of  the  photon  burst  is  included  to  enhance  the  resolution  of  the  image 
tube  and  that  it  would  be  mechanically  difficult  to  align  the  masks  and 
maintain  their  alignment.  In  response  to  the  first  criticism  it  shouid 
be  noted  that  the  binary  bar  masks  should  obviously  be  organized  in 
Gray  code.  The  periodicity  of  the  least  significant  bit  of  Gray  code 
happens  to  be  twice  that  of  ordinary  binary  so  that  the  photon  burst 
may  actually  be  two  or  even  three  resolvable  pixels  in  size  before  the 
bit  loses  its  significance.  In  response  to  the  second  criticism  the 
multiple  images  may  be  obtained  with  a  solid  glass  prism  on  which  the 
masks  may  be  firmly  glued,  thereby  assuring  stability  of  their  apparent 
relative  positions.  That  prism  also  leads  to  a  convenient  compact 
design  whose  section  is  shown  in  Figure  3. 

The  prism  is  fabricated  by  cutting  an  ordinary  right  angle  prism 
into  pie  shaped  slices.  These  slices  are  cemented  together  in  a 
circle  as  they  would  lie  in  a  complete  pie.  The  masks  are  glued  on 
as  indicated  in  Figure  A  with  the  result  that  the  virtual  images  of 
the  masks  are  superimposed.  The  effective  glass  thickness  beneath 
which  the  virtual  images  appear  would  ordinarily  contribute  substantial 
spherical  aberration,  but  that  aberration  is  suppressed  by  cementing 
a  plano-convec  lens  on  the  prism.  The  lens  is  chosen  so  that  the  virtual 
images  are  aplanatically  immersed  as  well  as  magnified  by  the  refractive 
index  of  the  glass.  The  resultant  prism  is  used  in  conjunction  with  a 
transfer  lens  that  gives  a  real  conjugate  image  of  the  superimposed 
virtual  images  of  the  masks.  The  phosphor  output  screen  of  the  image 
intensifier  is  to  be  located  at  that  real  conjugate  image.  Mask  align¬ 
ment  is  easily  accomplished  by  viewing  and  adjusting  the  registration  in 
the  real  conjugate  image  prior  to  the  setting  of  the  glue. 
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SPECIFICATIONS 

The  device  as  constructed  employs  15  standard  half-inch  diameter 
bialkali  photomultipliers  to  give  a  128  x  128  pixel  format.  Each  photo¬ 
multiplier  is  followed  by  a  one-transistor  emitter  follower  preamplifier 
and  a  high-speed  comparator.  The  comparators  each  have  an  adjustable 
reference  threshold  and  a  capacitor  feedback  to  stretch  the  pulses  to 
about  250  ns.  The  comparator  from  the  open  photomultiplier,  after  a 
delay  of  about  100  ns,  initiates  a  monostable  whose  leading  edge  serves 
as  clock  to  latch  the  outputs  of  all  the  other  comparators,  thus  registering 
coincidences  and  non-coincidences.  The  monostable,  whose  duration  is 
adjustable,  precludes  any  response  to  additional  pulses  during  that 
duration.  The  latched  outputs  are  then  available  for  digital  recording 
and  also  go  to  a  Gray  to  binary  analog  converter  composed  of  exclusive 
or  gates  and  thence  to  digital  analog  converters,  one  for  each  coordinate. 
The  D-A  converters  are  arranged  to  give  6-bit  magnitude  and  sign  so  that 
the  analog  outputs  are  bipolar  about  an  origin  at  the  center  of  the  picture 
format . 

The  image  tube  is  a  standard  2nd  generation  image  intensifier  except 
with  a  faster  phosphor.  The  model  we  have  is  a  Varo  3603  with  a  P22b 
rather  than  a  P20  phosphor,  and  its  photon  bursts  are  quite  readily 
detectable.  The  time  constant  of  the  phosphor  is  in  the  neighborhood 
of  1  microsecond,  even  though  the  decay  is  not  strictly  exponential. 

More  recently  we  had  the  opoortunity  to  borrow  a  Varo  3603  with  a  P46 
phosphor  whose  output  flashes  are  somewhat  shorter  with  less  of  an  after¬ 
tail  and  roughly  10  times  as  bright  as  those  of  the  P22b. 

The  remaining  problem  is  that  the  coincidence  detection  remains  un¬ 
reliable.  The  causes  of  the  unreliability  are  not  yet  identified;  they 
may  be  sensitivity  difficulties,  timing  difficulties  or  optical  difficul¬ 
ties.  As  an  aid  to  tracking  down  the  errors,  we  have  constructed  an  LED 
(green)  flasher  that  provides  300  ns  flashes  at  12  microsecond  intervals, 
and  the  brightness  of  the  flashes  is  step  selectable  over  the  range  to 
emulate  the  phosphor  flashes  from  the  intensifier  and  up  to  several  times 
brighter.  Of  course,  the  advantage  of  test  flashes  for  diagnosis  is  that 
the  location  and  timing  of  the  flashes  is  fixed. 
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COMMERCIALLY  AVAILABLE  DEVICE 

In  August  1979,  Surface  Science  Laboratories  of  Palo  Alto  announced 
an  imaging  electron  detector  for  sale.  Their  brochure  (enclosed,  in  part, 
as  Appendix  B)  indicated  that  a  photon  sensitive  version  was  under  develop¬ 
ment,  simply  by  placing  their  electron  detector  in  a  vacuum  envelope  behind 
a  photocathode.  Remembering  that  the  aim  of  our  project  is  speckle  imaging 
rather  than  detector  development,  and  also  in  view  of  the  unusual  time 
constraints  of  the  present  contract,  our  most  productive  course  has  been 
to  procure  the  commercial  detector.  In  essence  the  Surface  Science 
detector  is  a  carefully  tailored  resistive  anode  whereby  the  ratios  of 
discharge  distribution  among  four  electrodes  specify  in  analog  form  the 
two  coordinates  of  an  impinging  electron  burst.  In  a  sense  their  device 
is  the  electron  analog  of  the  simple  vignetting  photon  sensor  described 
earlier,  and  their  device  should  be  subject  to  similar  statistical 
limitations  of  precision.  They  accomplish  the  requisite  high  gain  by 
using  two  microchannel  plate  electron  multipliers  in  tandem.  At  any 
rate,  their  device  functions  quite  respectably  as  is  evidenced  by  their 
demonstrated  pictures. 

Development  of  the  Lockheed  photon  counter  will  nevertheless  continue 
because  the  direct  digital  approach  offers  greater  future  potential  in 
terms  of  resolution,  speed,  dynamic  range  and  format  versatility. 

SPECKLE  IMAGING:  AMPLITUDE  COMPENSATION 

Returning  to  the  topic  of  the  speckle  imaging  a  few  theoretical 

developments  have  arisen  since  publication  of  the  Applied  Optics  article 

(Appendix  A).  Dr.  T.  Tarbell  has  offered  a  compelling  argument  that  the 

amplitudes  of  the  spatial  frequency  components  are  necessarily  degraded 

as  a  result  of  the  image  spread.  Consider  the  gedanken  experiment  of 

mathematically  phase  shifting  the  spatial  frequency  components  of  a  point 

2 

image  until  the  point  spread  function  is  N  times  broader  (N  times  larger 
in  area).  According  to  Parseval's  theorem  the  integrated  power  of  the 
spatial  frequencies  would  remain  invariant  so  that  the  peak-to-peak 
speckle  amplitudes  would  diminish  by  N.  That  factor  is  inconsistent 
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with  conservation  of  energy  which  requires  that  the  envelope  peak  diminish 
2 

by  N  .  Consequently  the  atmospheric  spreading  must,  in  fact,  degrade  the 
amplitudes  of  the  spatial  frequency  components,  and  it  will,  therefore,  be 
necessary  to  calibrate  and  restore  the  degradation  according  to 

_  ^measured 

corrected  A  ...  , 

calibration 


where  A  calibration  is  the  measured  amplitude  for  the  spatial  frequency 
component  for  a  source  known  to  be  effectively  a  point.  As  is  well  known, 
when  the  degradation  is  severe  (small  ^caiibration^ '  t*ie  ^ar8e  required 
factor  of  restoration  mainly  amplifies  noise.  Wiener  filtering  provides 
a  form  of  compensation  that  takes  noise  into  account  so  as  to  eventually 
attenuate  rather  than  amplify  the  diminished  amplitudes.  Essentially, 
the  form  of  Wiener  filtering  is 


corrected 


A  A 

measured  calibration 

2  2 
A  +  o 

calibration 


where  a  is  a  constant  equal  to  the  rms  noise  present.  Asymptotically, 

when  A  .  >>a  then  the  naive  restoration  formula  is  followed, 

calxbration 

but  when  A  ,  .  «a  then  the  amplitudes  are  diminished  with 

calibration 

further  decrease  in  A  ,  .,  .  . 

calibration 

Now  the  speckle  imaging  algorithm  that  we  are  attempting  is  basically 
a  form  of  homomorphic  signal  processing,  so  that  we  are  concerned  and 
working  with  the  logarithms  of  the  spatial  frequency  components.  While 
it  would  be  awkward  and  time  consuming  to  reproduce  the  Wiener  filtering 
formula  in  logarithmic  form,  it  is  extremely  easy  to  compute  the  form 


loS  "corrected  '  log  "measured'  1050  '  T8  "calibration  '  l08C,|  ’ 

which  reproduces  the  asymptotic  behavior  of  the  Wiener  filtering  formula. 
Since  it  is  the  asymptotic  behavior  that  is  most  important,  this  much 

simplified  computing  procedure  should  be  perfectly  suitable. 
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SENSITIVITY  OF  SPECKLE  IMAGING 

There  are  also  further  thoughts  concerning  the  ultimate  sensitivity 
of  the  technique.  In  principle  the  reconstructed  picture  should  be 
reasonable  when  the  average  phase  settles  to  within  about  a  quarter  of 
a  turn  from  the  true  phase.  Suppose  that  photons  arrive  at  random 
positions.  Applying  the  algorithm  should  then  lead  to  a  random  walk 
of  the  phase,  and  a  random  walk  never  settles  to  a  legitimate  average. 
This  is  in  contrast  to  atmospheric  perturbations  where  shifts  must 
eventually  shift  back  to  present  a  proper  average.  At  the  other  extreme, 
when  there  are  lots  of  photons  representing  a  definite  picture  and  in 
the  absence  of  atmospheric  turbulence,  then  the  phase  measurements  will 
be  stable.  A  combination  of  the  extremes  such  that  the  random  perturba¬ 
tion  to  the  phase  never  leads  to  a  one-turn  slip  is  still  alright  for 
averaging.  But  if  the  perturbation  should  lead  to  a  slip  of  one  cycle 
in  the  middle  of  the  observing  run,  then  the  phase  average  will  end  up 
in  antiphase  with  the  desired  average.  In  other  words,  at  sufficiently 
low  brightness  levels  or  with  sufficient  sky  background,  the  photon  noise 
perturbations  can  destroy  the  validity  of  the  averaging  procedure.  These 
notions  relating  to  limiting  sensitivity  may  prove  to  be  tantamount  to 
just  another  point  of  view  on  calculations  that  others  have  already  made 
concerning  limiting  sensitivity  of  speckle  techniques.  The  additional 
point  of  view  is  still  welcome  if  it  can  simplify  our  understanding  of 
the  nature  of  the  limiting  sensitivity. 
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Introduction 

The  goal  of  this  work  is  to  attain  diffraction-limited 
performance  for  large  telescopes  in  spite  of  atmospheric 
turbulence.  That  possibility  has  become  realistic  as  the 
result  of  Labeyrie’s1  recognition  that  the  atmosphere 
rearranges  rather  than  obliterates  the  details  of  an 
image.  The  Abbe  perspective  that  views  an  image  not 
as  an  ensemble  ot  pixels  but  as  an  ensemble  of  spatial 
frequencies  is  heuristically  preferable.  Within  an  iso- 
planatic  field  (~3  sec  of  arc)  the  predominant  effect  o.‘ 
atmospheric  turbulence  is  to  shift  the  position  (phase) 
of  any  individual  spatial  frequency.  The  average  of  a 
bar  pattern  (spatial  frequency)  that  is  shifting  by  more 
than  half  a  bar  washes  out  the  contrast  of  that  pattern 
and  is  therefore  to  be  avoided.  The  essence  of  the  al¬ 
gorithm  to  be  described  is  to  average  separately  the 
strength  and  position  of  the  shifting  bar  patterns. 

Detection 

The  observational  raw  material  for  this  algorithm  is 
preferably  obtained  with  a  2-D  photon  counter,  which 
is  a  device  that  lists  both  coordinates  of  successively 
arriving  photons.  Such  a  list  is  clearly  more  informative 
than  an  image  because  the  list  may  be  irreversibly 
transformed  to  an  image.  Furthermore,  at  low  light 
levels  the  list  becomes  more  efficient  than  an  image 
because  at  low  light  levels  most  of  the  capacity  of  the 
image  is.  then  wasted  indicating  only  the  absence  of 
photons.  Given  the  motivation,  several  elegantly 
simple  schemes  for  constructing  such  a  device  haw  I  wen 
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devised  but  will  not  be  described  here  Photon  rates 
in  the  neighborhood  of  100- 10,000  photons/sec  in  the 
entire  picture  will  be  considered. 

Algorithm 

The  memory  si/e  of  the  computer  shall  he  four  times 
the  number  of  pixels  in  the  final  image.  Thus  a  I6A'  x 
16  bit  computer  will  suffice  in  this  example  to  yield 
64  X  64  picture  (4 K  pixels).  The  list  of  photon  coor 
dinates  is  off  line,  preferably  on  magnetic  tape.  The 
computer  memory  shall  be  divided  into  eight  equal 
sections  of  2 K  words  each.  The  first  section  will  include 
the  program,  which  is  very  short,  a  320-word  table  of  1  v.t 
cycles  of  a  sine  wave,  and  a  stack,  called  the  K  stack,  of 
length  ranging  from  2  to  256  words.  This  K  stack  will 
be  used  to  store  the  coordinates  of  photons  currentlv 
being  digested  by  the  computer.  The  length  of  the  h 
stack  shall  be  such  that  all  K  photons  arrive  while  the 
atmosphere  is  essentially  static,  i  e.,  within  about  1/50 
sec. 

Kach  remaining  section  is  to  be  construed  as  a  2  1) 
array  A,((',V),  with  (’  ranging  from  I  to  32  and  \ 
ranging  from  —31  to  +32.  These  arrays  will  be  initial 
i/.ed  to  zero.  1 '(imputations  will  then  commence  as 
follows. 

The  coordinates  (X  and  V  in  packed  form)  ot  the  first 
photon  are  read  in  to  the  first  word  of  the  K  stack. 
Copying  from  the  sine-cosine  tablp, 

.-I  j (/ V)  *  rust  V *  1  )  V  ). 

A  A l  V  i  =  ♦  >  V  i 

font  inning  to  read  in  photons  until  the  K  stack  is 

.ul, 

,1,iCVi  *  V  Y,r  *  Vi Vl  *  Rc|/, </  .Vi:  =  0,lC,V) 

K 

4  ..|  t  l  1  *=  -ini  V,  i  +  V.  Vi  --  tmi/,  1 1  *.  Vi;  -  h,  1 1  \  Vi 
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At  the  conclusion  of  this  first  filling  of  the  K  stack,?  = 

1 

In  this  fashion  wc  have  synthesized  the  Fourier 
transform,  /,(/  ',  VI,  of  the  i matte  corresponding  to  the 
first  K  photons  of  the  list,  and  a  few  comments  are  in 
order  The  t matte  is  necessarily  real,  and  therefore  the 
Fourier  transform  is  Hermit ian.  so  that  only  a  half  plane 
need  lie  stored.  In  other  words  the  combined  arrays  A  , 
and  A  •  map  reversibly  to  the  l>4  X  t>4  picture. 

Second,  it  K  is  small  this  synthesis  procedure  will  lie 
faster  than  a  fast  Fourier  transform.  Moreover,  the 
s|H-ed  is  still  greater  liecause  no  multiplications  at  all  are 
involved;  (A \l '  -I-  V,.  V)  is  obtained  from  additive  re¬ 
cursion.  sine  cosine  is  obtained  from  table  lookup,  and 
tile  coefficients  are  all  unity  liecause  each  photon  is  one 
unit. 

The  next  task  is  to  ascertain  the  strength  and  position 
ol  each  spatial  frequency.  Position  presents  the  awk¬ 
ward  problem  since  phase  as  determined  from  i/i  = 
arctaii  l  b/a\  is  only  a  partial  measure  of  position.  That 
phase  </>  is  only  the  fractional  part  of  a  turn,  integral 
turns  being  undetermined.  Dramatic  errors  will  occur 
if  we  attempt  to  average  a  set  of  values  knowing  only 
their  fractional  parts.  Therefore,  let  us  establish  a 
complete  position  that  includes  integral  turns, 

i  ii  i  i  i  i %  ,r: ::  +  a,  ,. 

with  =  ii.  The  bracketed  expression  is  reduced  to 
lying  between  ±!;.  turn  by  arithmetic  left  and  right 
shifts,  easy  computer  operations.  Basically  this  pro¬ 
cedure  tracks  the  phase  by  differentiating  and  inte¬ 
grating  so  as  to  follow  and  keep  track  of  integral 
turns 

Progression  in  time  occurs  as  follows.  The  K  stack 
is  full,  so  we  subtract  the  sine  and  cosine  terms  due  to 
the  initial  photon  from  A|  4f V’)  That  leaves  room 
in  the  K  stack  to  add  the  next  photon  so  that  we  now 
ha\eat?  =  2  t  he  Fourier  transform  of  photons  2  to  K  + 

I  mi  place  of  the  previous  Fourier  transform  of  photons 
I  to  A  It  is  this  stage  that  gives  rise  to  the  descriptive 
phrase:  photon  by  photon.  Actually  it  is  permissible 
to  shift  In  several  photons  as  long  as  several  is  appre 
i  mill v  less  than  K.  The  new  Fourier  transform  is  very 
rediind.ini  with  the  previous  one,  and  therefore  we  can 
be  coni ideiit  that  the  phases  will  have  shifted  by  less 
than  halt  a  cycle,  so  that  the  phase  t racking  is  legitimate. 
Furthermore,  I  his  recursive  progression  of  t  lie  Fourier 
t  raiisform  is  si  ill  taster  I  ban  starting  from  scratch  be 
cause  even  fewer  than  K  photons  are  involved. 

It  is  now  possible  to  consider  averaging  the  st  rength 
and  posit  ion  o|  t  he  components  so  ( hat  in  the  course  ol 
1  lie  i  oiiipnt.it n at  we  can  let 

.i,i r  i  I  s.  |/  1 1  . v  ,|. 

1  it  \  I  s_  t  it  i a 

where  |/|  In  A  h  '  =  a  cos +  b  sin,/,. 

At  the  i  oticlusion  ol  exhausting  the  catalog  7',  an 
estimate  ol  an  average  Fourier  transform  is  available 
as 


By  the  time  the  position  average  settles  to  within  %  turn 
accuracy,  this  estimate  may  be  reconstructed  back  to 
a  detailed  diffraction- limited  image  by  a  final  conven¬ 
tional  2  D  fast  Fourier  transform. 

Comparative  Comments 

At  this  stage  several  comparative  comments  are  in 
order.  The  first  is  a  comparison  with  Laheyrie’s  speckle 
interferometry.  Speckle  interferometry  works  for 
centrosvmmetric  objects.  In  other  words,  the  Fourier 
components  are  all  cosinusoidal  so  that  the  phases  may 
be  specified  as  zero  and  need  not  be  ascertained.  He 
uses  a  power  spectrum  averaging  for  the  strength  of  the 
components.  The  [lower  sped  rum  gives  slightly  dif¬ 
ferent  weighting  than  the  amplitude  average,  but  nev¬ 
ertheless  the  estimate  is  not  grossly  different. 

The  phase-tracking  procedure  may  be  also  compared 
with  the  Knox -Thompson-  algorithm.  The  Knox- 
Thompson  algorithm  differentiates  with  respect  to 
spatial  frequency,  rather  than  time,  for  tracking.  As 
a  result  their  procedure  can  o|ierate  with  only  occasional 
glimpses  of  the  scene.  Fewer  calculations  are  then  re¬ 
quired  under  relatively  high  brightness  conditions  be¬ 
cause  temporally  contiguous  glimpses  would  be  re¬ 
dundant.  The  Knox-Thompson  differentiation  pro¬ 
cedure  uses 

,  =  \/.,  |  |Z  ,|  .  ,i 

The  phase  difference  angle  is  always  less  than  Vj  turn 
so  that  the  imaginary  part  of  the  expression  is  an  am¬ 
plitude  weighted  approximation  to  that  difference 
angle.  Also  their  integration  procedure  accumulates 
round  off  errors.  The  tracking  difference  procedure 
may  he  applied  equally  in  spatial  frequency  with  only 
modest  extra  calculational  effort  by 

♦.  -  l*>.  A  !  | '  |  *  A.  ,.  A.  =11. 

A  summary  comparison  ol  the  estimated  phase  angles 
is 


where  =  <*.  and  t  is  just  a  dummy  index  for  the 
temporal  integration.  For  the  Knox-Thompson  algo¬ 
rithm 


and  ,  =  0.  When  the  phase  differences  arc  small,  the 
approximation  is  close,  hut  I  contend  that  the  expres¬ 
sion  (21  is  actually  more  legitimate  than  Kq.  (1).  The 
establishment  of  a  null  border  around  the  picture  would 
serve  to  diminish  the  phase  difference  angles  and  im¬ 
prove  the  approximation  of  Kq .  (1).  In  either  case 
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summation  with  respect  to  T  before  with  respect  to  V 
is  particularly  susceptible  to  accumulated  round-off 
error.  Interchanging  the  summation  leads  to  a  possible 
third  version, 

<*,  >  =  L  f  L  1<A,.,  ~  </><„- 1>]. 

/  t-»  Li  »i  J 

that  should  be  less  susceptible  to  round-off  even  though 
it  requires  more  arithmetic  operations. 

A  third  comment  relates  to  the  length  of  the  K  stack. 
As  stated  previously,  K  should  be  no  longer  than  the 
coherence  time  of  the  atmosphere,  in  other  words  less 
than  a  half  cycle  of  the  highest  really  significant  fluc¬ 
tuation  frequency.  That  time  depends  somewhat  on 
atmospheric  conditions  but  is  typically  about  1/50  sec. 
At  the  other  extreme  how  short  may  K  be?  The  answer 
is  2.  because  if  K  were  to  he  only  1,  the  reconstructed 
picture  would  consist  of  a  single  point  at  the  mean 
photon  arrival  position. 

A  fourth  comment  relates  to  the  susceptibility  of  all 
these  algorithms  to  telescope  aberrations,  particularly 
to  nonsymmetric  aberrations.  Consider  at  first  a 
nonturbulent  or  static  atmosphere.  Except  for  photon 
statistics,  all  glimpses  of  the  picture  are  then  identical, 
and  the  averaging  simply  reconstructs  that  same  pic¬ 
ture.  which  includes  telescope  aberrations.  If  we  now 
let  the  atmosphere  become  gradually  more  turbulent, 
more  variance  of  <l>  develops,  but  the  average  still  centers 
around  the  systematic  value  of  <t>  given  by  the  telescope 
aberrated  picture.  Therefore.  I  contend  that  these 
speckle  imaging  procedures  give  the  telescope  alierrated 
picture  rather  than  an  unadulterated  diffraction  limited 
picture.  On  the  other  hand,  it  should  not  he  difficult 
to  calibrate  the  aberrations  by  observing  an  unresolved 
star  and  subtracting  the  systematic  'f>  of  the  aberrations 
from  those  of  the  unknown  picture  to  obtain  an  una¬ 
berrated  diffraction-limited  picture.  A~{U,V )  is  re¬ 
served  for  storage  of  the  reference  4>ref.  Such  calibra¬ 
tion  requires  only  that  the  telescope  aberrations  be  re¬ 
producible. 


Improvements 

Let  us  now  return  toward  improving  the  estimate  of 
the  algorithm.  The  atmospheric  effect  of  shifting  the 
complex  Fourier  components  is  mathematically  mul¬ 
tiplicative.  There  is  a  less  severe,  but  nevertheless 
appreciable,  effect  on  the  Fourier  amplitudes  that  also 
is  a  multiplicative  process  and  must  lie  averaged  out.  It 
should  be  preferable,  however,  to  have  an  additive 
process  so  that  the  fluctuating  component  may  have  a 
proper  zero  mean.  The  requisite  process  conversion  is 
accomplished  simply  by  taking  logarithms  of  the  com¬ 
plex  Fourier  components,  a  procedure  that  is  used  i  i 
cepstral  analysis'  for  use  with  homomorphic  signal 
processing,  of  which  this  speckle  imaging  may  he  a  2-1) 
example.  If  V.  is  complex,  inZ  is  complex; 

/  ~  a  +  ih  -  tMpilrt/)  t-xpui  +  /,/). 

We  note  that  •>  =  ln|Z|  and  ti  =  2ir'l>,  so  that  it  should 
he  preferable  to  calculate 


A«(l M')  =  ln|4  1 


and  use  the  estimate 


/  |  /  v  | 

Est<Z)  =  exp  y-  V  In/,  j  =  exp  i/t,  +  2mA., I 


A  further  improvement  may  be  available  by  noting 
that  if  the  amplitude  of  a  Fourier  component  happens 
to  be  strong,  there  will  be  less  uncertainty  in  the  de 
termination  of  its  position.  Thus  an  amplitude 
weighting  is  desirable  for  the  position  average,  and  not 
to  discriminate  against  the  real  part,  that  also  should 
be  weighted.  Thus  we  should  calculate 


AtH’.V)  =  L  |4|  In  1 4 1. 


AM'.V)  =  {_  |4|*,. 


=  v  |4 


and  use  the  estimate 


Est(Z)  =  pxp 


v|Z|  \nZ 

i_|/f 


I 

exp  —  (.4*  +  2»i.4;,l. 


as  the  basis  for  reconstructing  the  picture.  Note  that 
the  entropy  form  of  the  expression  has  entered  in  a  very 
natural  way.  Insofar  as  entropy  is  identified  with 
probabilistic  information,  I  suspect  that  the  picture 
based  on  this  last  estimate  is  the  most  probable  picture, 
in  a  Bayesian  sense,  that  would  have  given  rise  to  the 
observed  photon  catalog. 


Speed 

For  a  minicomputer  using  software  multiply  the 
conversion  from  rectangular  (a,/))  to  polar  (|Z|.</,i 
notation,  and  the  logarithmic  function  takes  much  more 
time  than  the  synthesis  of  the  Fourier  transforms.  For 
the  64  X  64  scale  of  the  example  the  program  runs  about 
500  times  slower  than  real  time,  so  that  10  see  of  data 
would  require  about  1 h  of  computation.  That  rate 
is  not  at  all  discouraging,  and  much  faster  rates  (even 
real  time  for  TV  size  pictures)  may  he  anticipated 
eventually  because  the  algorithm  is  so  suitable  for  mixed 
parallel-serial  processing. 

The  attainment  of  diffraction-limited  resolving  power 
in  the  Rayleigh  sense  means  that  the  spatial  frequencies 
are  each  correctly  positioned  to  within  about  V4  bar 
(fringe).  Thus  if  the  rms  atmospheric  phase  pertur¬ 
bation  is  five  fringes  we  may  expect  to  need  about  20- 
=  400  statistically  independent  samples  for  their  aver 
age  to  settle  within  that  Rayleigh  criterion.  Presuming 
that  samples  separated  hv  about  1/10  see  are  fairly  in¬ 
dependent,  about  40  sec  of  observation  should  he  suf¬ 
ficient  to  quell  the  effects  of  atmospheric  turbulence. 


Noise 

Let  us  now  examine  the  behavior  of  speckle  imaging 
with  respect  to  noise.  One  principal  task  of  a  large 
telescope  is  to  seek  ever  fainter  stars.  The  obstacle  is 
the  background  illumination  from  the  night  sky,  so  that 
the  function  of  the  telescope  is  not  only  to  gather  pho- 
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tons  but  also  to  isolate  the  stellar  photons  as  thoroughly 
as  possible  from  the  night  sky  photons.  Superior  iso¬ 
lation  is  in  fact  one  major  goal  of  the  space  telescope. 
For  speckle  imaging  a  uniform  background,  along  with 
its  statistical  fluctuations,  gets  rearranged  into  a  uni¬ 
form  background  with  comparable  statistical  fluctua¬ 
tions.  A  star,  on  the  other  hand,  gets  rearranged  from 
an  extended  patch  to  a  concentrated  point.  Therefore, 
the  contrast  of  the  star  to  background  should  be  aug¬ 
mented  highly,  accomplishing  from  the  ground  that 
major  goal  of  the  space  telescope. 

The  repercussions  with  respect  to  very  large  tele¬ 
scopes  should  also  be  important.  The  motivation 
toward  larger  apertures  should  now  be  enhanced  by  the 
new  opportunity  to  increase  resolving  power  as  well  as 
collecting  area.  Economy  dictates  that  very  large  ap¬ 
ertures  be  assembled  from  pieces  with  the  result  that 
random  aberrations  will  be  much  larger  than  the  dif¬ 
fraction  limit  even  though  the  quality  commensurate 
with  atmospheric  seeing  is  attained.  However,  those 
aberrations  must  be  reproducible  if  they  are  to  be  re¬ 
moved  by  calibration.  The  motions  of  any  assembly  of 
servo-controlled  pieces  will  certainly  hinder  greatly,  and 
probably  rule  out  completely,  this  use  of  speckle  imag¬ 
ing  to  gain  the  diffraction  limit.  It  is  far  more  likely 
that  an  immobile  aperture,  such  as  the  Arecibo  concept 
that  1  have  long  advocated,4  might  succeed.  A  very 
narrowband  (~1-A)  color  filter  would  permit  operation 
independent  of  primitive  phasing  (piston)  errors  of  the 


pieces.  It  should  be  remembered  that  the  decrease  in 
flux  by  the  restriction  to  narrow  bandwidth  also  applies 
to  the  night  sky  and  furthermore  would  be  compensated 
by  the  enormous  collecting  area  of  the  aperture. 

Remarks 

It  is  a  pleasure  to  call  attention  to  a  commendably 
clear  and  thorough  review  article  by  Worden6  describing 
previous  work  in  speckle  imaging.  The  present  work 
originated  while  1  was  at  the  Smithsonian  Astrophysical 
Observatory;  ordinarily  publication  should  have  been 
deferred  pending  experimental  verification.  However, 
the  necessary  experiments  have  had  to  be  adjourned 
since  I  am  no  longer  employed  there.  Nevertheless  the 
goal  of  attaining  diffraction-limited  performance  for 
large  telescopes  seems  so  worthwhile  and  so  accessible 
that  1  still  hope  for  the  opportunity  to  resurrect  the 
experimental  work. 

Postscript.  Thanks  to  my  new  employer,  this 
speckle  project  is  being  started  afresh. 
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APPENDIX  B 


JUST  ONE 

electron/msec 

...HONEST. 

That’s  the  Illumination  we  used 
to  take  this  picture 

With  our  unique 

IMAGING  PARTICLE  DETECTOR 


Each  dot  on  this  CRT  photo  represents  a  single 
electron  passing  through  an  object  mask  onto 
our  detector  (30  sec  exposure) 


USES  •  Detector  for  electron  specttometers;  mass  spectrometers; 

UV/X-Ray  imaging  systems  particle  physics  experiments  etc 

FEATURES  •  High  total  quantum  efficiency 

•  Useful  in  fluxes  up  to  105  events/sec  dead  time  7Msec 

•  Excellent  spatial  resolution  >100  line  pairs/in 

•  Compact:  2"  Dia  .,  1.5"  thick  including  electron  multipliers 

•  Extremely  rugged,  stable,  and  long-lived 

•  Digital  or  analog  outputs  available 

HOW  IT  WORKS 


j-h 


Position 

computing 

electronics 


X 

Y 


incident  flux 
(particles,  photons) 

electron  multipliers 
(gain  m  10') 

resistive 

anode 

to  CRT, 
computer,  etc 


Position  signals 


•  A  part'cle  strikes  the  electron  multipliers  producing 
a  pulse  of  106  electrons 

•  This  charge  is  collected  by  4  electrodes  on  a 
resistive  anode 

•  The  x,y  position  ot  the  particle  is  computed  from 
the  charge  division  among  the  electrodes 


Call  Mike  Kelly  at 

•  URPACi  SCIINCE  LABORATORIES 

4151  Middleheld  Road 
Palo  Alto,  CA  94303 
415  493-0229 


MicroChannel  Plate 
A  ssembly 
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APPENDIX  B  (cont'd) 


PERFORMANCE 


The  two  photographs  on  the  next  page  show  the  performance  of 
the  overall  two-dimensional  system.  They  are  photos  of  a  storage 
scope  connected  to  the  "X"  and  "Y"  outputs  of  the  electronic  unit. 

The  input  of  the  two  microchannel  plates  was  covered  by  a  mask  having 
the  test  pattern  on  it  and  illuminated  by  a  source  of  electrons. 
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APPENDIX  R  (cont'd) 


Performance  of  Two-Dimensional  Resistive  Anode  Detector 


Figure  1.  3  Example  of  detector  linearity:  electron  image  of  shadow  mask 
with  0.  1  inch  hole  spacing. 


Figure  1.4  Example  of  detector  resolution:  image  size  is  0.  1  inch/division 


APPENDIX  B  (cont'd) 


Optical  Detectors 


To  use  the  imaging  detector  for  visible  radiation,  a  photo- 
cathode  must  be  added  as  shown  in  the  figure  below.  SSL  does  not 
at  present  have  such  a  configuration  available,  but  one  is  under¬ 
development  . 

Visible/UV  radiation 


For  further  details  please  contact  Mike  Kelly  at  Surface 
Science  Laboratories 
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